Quantum entanglement in dense multiqubit systems 
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The pairwise entanglement of an arbitrary atomic pair randomly extracted from a laser-driven 
dense multiqubit sample in the presence of quantum dissipation due to spontaneous emission is 
considered. The dipole-dipole interaction between the particles shifts the laser-qubit resonance 
frequency and consequently modifies the quantum entanglement. By means of an appropriate tuning 
of the laser frequency, one can optimize the entanglement in this system. For large ensembles, the 
maximum entanglement occurs near the laser parameters where the steady-state sample exhibits 
phase transition phenomena. 
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I. INTRODUCTION 

The existence of entanglement in multiparticle sam- 
ples is an important fundamental problem However, 
this topic is difficult due to the fact that a general cri- 
terion that quantifies the entanglement in such systems 
does not exist. One may consider the state that can- 
not be represented as a factorized product of individual 
states describing each particle separately as an entangled 
multiparticle state. As a consequence, the condition for 
entanglement of formation of an arbitrary state of two- 
qubits and an inseparability criterion for continuous 
variable systems Q were proposed. Other criteria are, 
respectively, the entanglement of distillation Q, the rel- 
ative entropy of entanglement [5[ or the so called nega- 
tivity |6(. Based on these formulations, numerous inves- 
tigations were performed in order to describe the entan- 
glement in various systems. In particular, entanglement 
induced by a single-mode heat environment was shown 
to occur in Q while entanglement between two atoms 
in an overdamped cavity injected with squeezed vacuum 
can be effectively created [8] . Two initially entangled and 
afterward not interacting qubits can become completely 
disentangled in a finite time 0]. However, dark periods 
and revivals of entanglement in a two-qubit system was 
shown to occur via environmental vacuum modes [lol ]. 
An experimental demonstration of continuous variable 
entanglement using cold atoms in a high-finesse optical 
cavity was presented in while experimental evidence 
of quantum entanglement of a large number of photons 
was proved in Ref. [l^ . Also, the entanglement of for- 
mation for an arbitrary two-mode Gaussian state was 
demonstrated in jl3j |. 

Usually, many-particle systems show critical behaviors 
such as quantum phase transitions [lil . [l5l . The en- 
tanglement properties of the quantum phase transition in 
particular systems [l7l [l8j as well as their instability [l9j 



have been investigated in detail. There, the concurrence 
[2], which characterizes the entanglement between two 
spins (i.e. pairwise entanglement) after tracing out other 
all others, was calculated. This quantity is well-suited for 
multiparticle systems since it does not depend on the two 
spins selected because all spins are completely equivalent 
[20(. Moreover, arbitrary symmetric multiqubit states 
are s pin squeezed if and only if they are pairwise entan- 
gled [21(. Note that spin squeezing was experimentally 
measured in [H, [23[ . 

Therefore, in this paper we investigate the pairwise en- 
tanglement in a laser-pumped dense multiparticle sam- 
ple which may characterize the quantum nature of the 
system as a whole. The collectivity among the par- 
ticles is mediated through the surrounding electromag- 
netic field vacuum reservoir [24j |. This system is known 
to exhibit first- or second-order phase transitions in the 
steady-state. We estimate the concurrence for an arbi- 
trary atomic pair randomly extracted from the atomic 
ensemble. The influence of the dipole-dipole interaction 
among particles on the entanglement creation is discussed 
in detail. In particular, the dipole-dipole interactions sig- 
nificantly enhance the pairwise entanglement in a two- 
qubit system. For larger samples, the maximum possible 
value of the concurrence is less sensitive to the dipole- 
dipole interactions. The relationship between the entan- 
glement and phase transition phenomena in this system 
is established. We found that the maximum of the pair- 
wise entanglement occurs near the critical steady-state 
behaviors where phase transitions take place for larger 
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FIG. 1: Schematic of a laser-pumped multiqubit two-level 
ensemble. Our interest is whether or not an arbitrary pair is 
entangled inside the sample. 
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samples. 

The paper is organized as follows. In Section II, we in- 
troduce the system of interest and solve it in the steady- 
state. Section III deals with the pairwisc entanglement in 
a symmetric multiparticle sample. Section IV describes 
the obtained results. We finalize the article with conclu- 
sions presented in Section V. 



II. SYSTEM OF INTEREST 

We consider a dense multiqubit two-state (|e) and \g)) 
sample pumped with a coherent laser field and damped 
via its interaction with the environmental electromag- 
netic field reservoir (see Fig. [T]). The frequencies of all 
atoms are identical and equal to ujq, while the laser fre- 
quency is u>l- The particles are placed in a volume with 
dimensions small compared to the emission wavelength. 
As a consequence, we assume that the dipole-dipole in- 
teraction potential, 5, is identical for all qubit pairs. 
App ropriate systems could be driven solid-state samples 
[25j ] . superconducting qubits pumped multiparti- 

cle systems trapped inside optical cavities [28] or double 
quantum dot excitonic systems 29]. Such a model can 
be described in the Born-Markov, the rotating-wave and 
the electric-dipole approximations by the master equa- 
tion approach [30J: 



dt P 



i[H ,p] - j{[S + ,S-p] + [ P S + , S-}}, (1) 



where 



Hq = AS Z - 5S + S~ + n(S + + S~), 



(2) 



with A — u>q — u>l + 8 = A + 5. Here, f2 is the Rabi fre- 
quency and 7 is the single-particle spontaneous emission 
decay rate. The collective angular momentum operators, 
S' ± and S z , are defined in the following way: 



5 + 



N 



S~ 



(S 



+\1 



N 



(3) 



where the raising and lowering operators for each qubit 
are denoted by cr+ = \e)jj(g\ and aj — (<7y~) while 
a Z j = \e)jj(e\ — \g)jj(g\- The collective operators 
obey the commutation relations for su(2) algebra, i.e. 
[S z , S ,± ] = ±S ± and [S+, S~] = 2S Z . 

Since our main interest here is the steady-state entan- 
glement for random pairs of qubits extracted from the 
whole ensemble, we present the steady-state solution of 
the master equation ([T]) (3l| : 



JY 



= z- x £ C nm (S-) n (S+) 



n.m— 



where C nm = cr n \a*Y m a nm with 

_ r(l +n + P)T(l +m + (3*) 
anm ~ n!m!r(l + /3)r(l + /3*) " 



(4) 



(5) 



Further, a = ifl/ (7 + iS) and f3 = iA/(7 + iS) while Z is 
chosen in such a way that Tr{p s } = 1. The trace can be 
performed using the relations 



S+\S, I) = y/(s-l)( 8 + l+l)\8, l + l), 
S-\S,1) = y/(s + l)(3-l+l)\s,l- 1), 



(6) 



where the collective Dicke states \s, I), with s = N/2 and 
—s<l<s, are the eigenstates for the operator S z and 
the operator of the total "spin" S 2 : 



S z \s,l) = l\s,l), 
S 2 \s 7 l) = s(s + l)\s,l). 



(7) 



Therefore, 



v - ,_ 2n (iV + n+l)!(n!) 2 

Z = 2_^a nn \a\ — — mM/0m | — r (8) 



n=0 



(N — n)\(2n + 1)! ' 



The steady-state atomic variables of interested can be 
obtained using Eqs. (J3HHJ) : 



N 



((S+T(S z nS-y) = Z- 1 C n-f,n- P 

n=max{f.p} 

N-n 



E 



(N — m)!(m + n)\ 
(N — m — n)\m\ 



(N/2-m-n) r , (9) 



where {p, r, /} G {0, 1, 2, — , N}. The exact steady-state 
results for large samples show that, for A = S = 0, 
the steady state averages are a continuous function of 
the driving field parameter 2fl/(N^) but their deriva- 
tives with respect to that parameter are discontinuous at 
2Vl/{N r y) = 1. This behavior is reminiscent of a second- 
order phase transition. If A ^ but 6 = there is 
not such a critical phenomenon (l5j . However, when 
{A, 6} 7^ the steady state solution ^ predicts a first- 
order phase transition for large systems (l5l . |3~0 ] . 

Further, we shall use the expectation values of the col- 
lective operators given by Eq. © to estimate the pair- 
wise entanglement among any pairs randomly extracted 
from the multiqubit ensemble. The relationship between 
the entanglement and above mentioned phase transition 
phenomena will be established. 



III. PAIRWISE ENTANGLEMENT 

For a mixed state of qubits {a, b} e {1, • • ■ , N} with 
density matrix p a b, the concurrence C is defined as 



C = max{0, Ai — A2 — A3 — A4}, 



(10) 



where the quantities Ai (i <E {1,2,3,4}) are the square 
roots of the eigenvalues of the matrix product 



R = Pab(cr ay ® O h y)p* ah {o ay ®O by ), 



(11) 
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FIG. 2: (color online) The concurrence C as function of 
2Q/(N~i) for N = 2 and A + S = 0. The solid curve is 
for 5/7 = 0, the long-dashed line corresponds to 5/j = 5, the 
dashed-dotted line is for S/j = 10 while the dotted curve is 
for 5/7 = 15. Thin solid lines show corresponding results for 

C^Jf defined in Eq. ||I7). 



in descending order. Here, p* h denotes complex conjuga- 
tion of pab, and Pauli matrices for the two-level 
systems (i £ {a, b}j. The values of the concurrence range 
from zero for an unent angled state to unity for a maxi- 
mally entangled two-particle state Q. Note that we do 
not rescale the obtained concurrence with the number of 
atoms N, in contrast to some previous works (l8l. [20|. 

The state of two qubits which is symmetric under the 
exchange of the sub-systems can be represented in the 
basis {\e a e b ), \e a g b ), \g a e b ), \g a gb)} as follows: 



Pab 



( Pll Pl2 Pl3 Pli 

pIX P22 P23 P2i 

P31 P32 P33 P3A 

V Pil PA2 PA3 P44 



(12) 



For symmetric multiparticle systems like ours, the den- 
sity matrix p a t and hence the concurrence do not depend 
on the specific choice of qubits a and b. Therefore, the 
matrix elements of p a b can be determined by: 



Pn 



N 2 — 2N + 4(S 2 ) + 4(7V - 1){S Z ) 



Pl2 = Pl3 



4N(N - 1) 
N(S+) +2(S+S Z 



2N(N- 1) 



Pu 



N(N — 1) 
N 2 



P22 — P23 



P21 
MS 2 ) 



(P12)*, 



P24 

P31 
P34 
P43 

PU = 



AN(N — 1) ' 
(S+)(N -2)-2(S+S z ) 
2N(N-1) 

(Pia)*, P32 = (P23 
P24, P41 = {puT, 

(Pm)*, 

N 2 — 2N + 4:(S 2 ) -A(N 



P33 — P22, 
P42 = (p2i)*, 



AN(N-l) 



(13) 




FIG. 3: (color online) The concurrence C as function of Q/j 
and A/7. Here N = 2 and 5/7=5. 



As it can be observed here, the two-particle density ma- 
trix is expressed in terms of expectation values of col- 
lective operators. The approach is valid when one has a 
symmetric state of N two-level particles, where by sym- 
metric, we mean symmetry under any permutation of the 
qubits [13, HH- Thus, the relations pfB|) will help us to 
estimate the pairwise entanglement of an arbitrary qubit 
pair inside the multiparticle ensemble. 

We also note that it was found in the literature that 
for density matrices which due to the studied model are 
restricted to the simpler case P12 — P13 — P24 — P34 — 0, 
a compact analytic expressions for the concurrence can 
be found, which is given by [32| 

a e /=max{0,cW,C^} , (14) 
C« =2(\p 14 \-^psT), (15) 
a^ = 2(|p 23 |-VpiIpiI). (16) 

We will see that which for the density matrix ele- 
ments in Eq. (|T3|) evaluates to 



C 



(i) 

re/ 



2(|pm| - P22) 



(17) 



will turn out to fit most of our results well. In the next 
Section, we turn to a characterization of the two-qubit 
concurrence in various parameter ranges. 



IV. RESULTS AND DISCUSSION 

In Figure El we plot the steady-state concurrence for a 
driven two-qubit system. The solid curve describes the 
case when dipole-dipole interaction among the particles is 
absent while the laser field is on resonance with the qubit 
frequency. The other curves characterize the situations 
where the dipole-dipole interaction increases while the 
laser field frequency obeys the condition A+<5 = 0. These 
behaviors can be understood in the following way. In the 
limit of small interparticle distance, the two-particle sys- 
tem is equivalent to a three-state ladder system: \g a gb) 



4 



0.06 F 




FIG. 4: (color online) The pairwise entanglement C as func- 
tion of 2Q/(N'y) for TV = 6. The solid red line is for 
8 = A = 0, the long-dashed black line stands for 28/ (TV7) — 1 
and 2A/(A r 7) = —0.25, the dashed-dotted blue curve for 
28 /(Nj) = 1 and 2A/(A?7) = -1.4, while the dotted green 
curve corresponds to 28 /(JV7) = 1 and 2A/(7V7) = —2. 
Thin solid lines show corresponding results for C^l defined 
in Eq. (fTT|) . 



0.012 F 




FIG. 5: (color online) The concurrence C as function of 
2fi/(/V7) for TV = 50. The solid red line is for 8 = A = 
0, the long-dashed black line stands for 25/(N r y) = 0.1 
and 2A/(TV7) = —0.1, the dashed-dotted blue curve for 
2<V(TV7) = 0.2 and 2A/(N'y) = -0.2, while the dotted green 
curve corresponds to 2S/(N-y) = 0.3 and 2A/(N'y) = -0.3. 
Thin solid lines show corresponding results for C^l defined 
in Eq. (TIT)) . 



being the ground state, |s) = (|e a <7{,) + \ebg a ))/y/2 is the 
intermediate level and the excited state is \e a ef,). The in- 
termediate level |s) is shifted by the dipole-dipole inter- 
action potential from the qubit frequency ujq. Therefore, 
in order to be on resonance A + S = 0, one has to modify 
the laser frequency as shown in Fig. [5] by the interrupted 
curves. As a result, the maximum of the concurrence 
moves to higher laser strengths. To further characterize 
the two-qubit system, in Fig. [31 we fix the dipole-dipole 
interaction potential, i.e. J/7 = 5, and change the fre- 
quency of the laser held. Interestingly, the concurrence 
can be further increased in comparison to the one shown 
in Fig.[^]by the long-dashed curve, where C max ~ 0.34 is 
achieved. For example, in Fig. [3] the maximum concur- 
rence is C max w 0.4 when A/7 = —12 and fi/7 m 1.8. 
Moreover, the concurrence shows oscillatory behaviors 
for larger field intensities as a function of A (see Fig. [3]). 

Further, we focus on the pairwise entanglement of a 
pair of qubits arbitrarily extracted from a multiparti- 
cle sample. We start with a moderate ensemble con- 
taining several particles. Figure (0| depicts the pair- 
wise entanglement for such a system with TV = 6 qubits. 
There, the solid line stands for 6 — A = while all 
others are for a fixed dipole-dipole interaction and vari- 
ous laser field detunings. Similar to the two-qubit case, 
the concurrence can be enhanced for moderate systems 
by modifying the external controllable laser parameters, 
i.e. {f2,A}, as well as the dipole-dipole interactions 6. 
For a fixed dipole-dipole interaction, the concurrence in- 
creases while changing the qubit-laser detuning and then 
decreases with an additional increase of the detuning (see 
Fig. |4|). For larger samples, i.e. TV ^> 1, the maxi- 
mum possible value of the concurrence is less sensitive 
to the dipole-dipole interparticle interactions, i.e. has al- 
most the same value as for the case of 6 = A = (see 



Fig. [5] and Fig. 0. Here, an important issue is that for 
S — A = and larger samples the maximum of the con- 
currence occurs near 2Q/(Nj) = 1 [TtJ and drops there 
abruptly (see Fig. [5] and Fig. [5]). At this critical point, 
i.e. 2n/(N"f) — 1 and TV ^> 1, there exists a second- 
order phase transition behavior of the collective atomic 
variables [l5| . In the presence of the dipole-dipole inter- 
actions, i.e. 5 =/= 0, the maximum of the pairwise entan- 
glement occurs at particular laser frequencies and may 
shift towards larger values of 2f2/ (TV7) depending on the 
laser-qubit detuning. Here again, the maximum of con- 
currence suddenly goes to zero at some critical values of 
2fi/(TV7) (see Fig. [5] and Fig. [5]). However, at these crit- 
ical values of the pumping parameter, a first-order phase 
transition of the collective atomic variables takes place 
[3l|. Therefore, one can conjecture that always for larger 
samples and in the presence of quantum dissipation due 
to spontaneous emission, the maximum of pairwise en- 
tanglement occurs nearby the critical behaviors of the 
collective atomic variables, i.e. phase transitions. 

It is interesting to see that in particular for larger sam- 
ples (TV > 1) and for parameters close to those which 
maximize the obtained concurrence, the reference 
in Eq. (jTTJ) fits our results very well. In Figs. [5] and [51 
no difference to the exact results can be seen, while the 
coincidence is equally good in Fig. [4] except for the dot- 
ted curve with detuning away from the value leading to 
the optimum concurrence. For the two-atom case, how- 
ever, C^~\ only reflects our numerical results for the case 

of absent dipole-dipole interaction, see Fig. [21 Note that 
(2) 

C^e/ does not indicate entanglement in any of our exam- 
ples. The good agreement with C^X suggests that it is 
in fact entanglement between the collective ground and 
excited states |4) and |1) which leads to the observed non- 
zero concurrence, rather than entanglement between the 
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FIG. 6: (color online) The concurrence C as function of 
20. /(N^) for N = 74. The solid red line is for 5 = A = 
0, the long-dashed black line stands for 28/{Nf) = 0.1 
and 2A/(A?7) = —0.1, the dashed-dotted blue curve for 
2S/(Nj) = 0.15 and 2A/(Nj) = -0.15, while the dotted 
green curve corresponds to 2S/(N r f) = 0.2 and 2A/(N / y) = 
—0.2. Thin solid lines show corresponding results for C^Q 
defined in Eq. JT7J ■ 

two intermediate states \2) and |3). This is somewhat re- 
markable, as the entangled symmetric state (|2) + |3))/v / 2 
naturally arises as an eigenstate of the coherent dynamics 
of dipolc-dipole interacting two-level systems. 



V. SUMMARY 



In summary, we investigated the pairwise entangle- 
ment of a randomly extracted qubit-pair from a laser- 
driven dense multiparticle ensemble. The whole system is 
damped collectively via its interaction with the environ- 
mental electromagnetic field vacuum modes. We started 
with a two-qubit system where the dipole-dipole inter- 
action significantly enhances the concurrence. Larger 
values for the concurrence in such systems can be ob- 
tained by adjusting the laser field detuning A around 
the frequency shift due to dipole-dipole interaction. In- 
tense driving laser fields lead to an oscillatory behavior 
of the concurrence in an off resonant two-particle sample 
as a function of the detuning. For moderate sample sizes, 
the concurrence can be further increased by suitably ad- 
justing the laser frequency and dipole-dipole interactions, 
similar to the two-qubit systems. Interestingly, the pair- 
wise entanglement goes abruptly to zero for larger sam- 
ples, at some particular values of the pumping parameter 
2fi/(A^7). At these critical values, phase transition takes 
place. 
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